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Abstract. The modelling of interface migration and the associated diffusion mechanisms at the nanoscale 
level is a challenging issue. For many technological applications ranging from nanoelectronic devices to 
solar cells, more knowledge of the mechanisms governing the migration of the silicon amorphous/crystalline 
O l interface and dopant diffusion during solid phase epitaxy is needed. In this work, silicon recrystallisation in 

the framework of solid phase epitaxy and the influence on orientation effects have been investigated at the 
atomic level using empirical molecular dynamics simulations. The morphology and the migration process 
of the interface has been observed to be highly dependent on the original inter-facial atomic structure. The 
[100] interface migration is a quasi-planar ideal process whereas the cases [110] and [111] are much more 
O , complex with a more diffuse interface. For [110], the interface migration corresponds to the formation and 

C/} ■ dissolution of nanofacets whereas for [111] a defective based bilayer reordering is the dominant re-growth 

process. The study of the interface velocity migration in the ideal case of defect free re- growth reveals no 
^ , difference between [100] and [110] and a decrease by a mean factor of 1.43 for the case [111]. Finally, the 

influence of boron atoms in the amorphous part on the interface migration velocity is also investigated in 
f ' the case of [100] orientation. 

PACS. 68. 55. A- Nucleation and growth - 81.15.Aa Theory and models of film growth 

Q\ 1 Introduction hand, the concept of a diffuse interface where an ordering 

Q\ . flux is associated to the interface migration is also pro- 

Advanced technological research work dedicated to silicon P osed - This flux is assumed to be generated by the coun- 

^ materials process engineering needs to control the silicon terbalance between the atomic forces interaction induced 

^sD recrystallisation process at the nanometer scale for nano- b y the short-range order and the thermal motion ©. In 

O electronic [H2] or photo- voltaic devices [3 4 . The silicon this framework, solid phase regrowth is often described 

^ film re-growth occurs during annealing when amorphous as an thermal activated process shown by the following 

T 7^ regions or layers are created by parasitic phenomena like in equation : 

^ . vapor-liquid-solid (VLS) grown silicon nanowires or inten- ^ 

' tionally, to limit diffusion effects like in advanced junction vspe(T) = vo exp [ -] (1) 

r/N depth processing [5 j. However, it is really difficult to reach 

q$ a deeper understanding of the many effects taking place where T is the interface temperature, vo is the ex- 
' during regrowth. For example, the rich interplay between ponential prefactor and E a is the activation energy. For 
the nucleation rate, the interface roughening, the defects silicon phase epitaxy, the debate focused on the nature 
generation or dissolution and the dopant diffusion is diffi- of defects (dangling bonds, kink sites, vacancies, inter- 
cult to characterise by using only experimental fabrication stitials) on their location (crystal, interface, amorphous 
and characterisation techniques. part) and on their influence on the limited growth rate 
From the theoretical point of view, there is a debate [8 . The influence of interface orientation, stress or dopant 
on the general physical description of crystal growth [6 . is also under debate. Atomistic process modelling [9] in 
The concept of a localised interface approach where the the framework of molecular dynamics is expected to help 
nucleation and growth process take place was first intro- to clarify the physical picture. This method could proba- 
duced by Burton, Cabrera and Frank [7 . On the other bly provide more insight into the re-growth mechanisms 

taking place at the amorphous/crystalline interface dur- 
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sociated diffusion effects during re-growth are much more 
complex to simulate compared to the bulk case [10] . More- 
over, information technology achieved significant progress 
and enabled the simulation of relatively large atomic con- 
figurations that could match experimental systems. 

Preliminary simulation studies of the silicon solid phase 
epitaxy related to junction processing by molecular dy- 
namics were performed by Mootoka et at. [11]. But many 
open questions remain on solid phase epitaxy like the sub- 
strate orientation influence despite clear technological in- 
terest. More knowledge is available in the case of melt- 
interfaces re-growth [12]. This problem has also been in- 
vestigated using a finite elements based method or lattice 
Kinetic Monte Carlo (KMC) [13] but numerous challenges 
are faced to take into account the complex reordering pro- 
cesses at the atomic level. The main objective of this pa- 
per is to present the work performed in order to investi- 
gate the influence of orientation and dopant effects on the 
amorphous silicon re-growth mechanisms. The first sec- 
tion briefly describes the computational details and the 
method used in this work. Some complementary tests re- 
lated to the influence of the initial properties of the dis- 
ordered layer are also provided. In a second part, most 
important results concerning the influence of the inter- 
face orientation on recrystallisation mechanisms are pre- 
sented and discussed. Preliminary results on the influence 
of Boron atoms on the interface migration velocity are also 
provided. Finally, some important conclusions are drawn 
in the last section. 



2 Methods 

2.1 Computational details 

Empirical molecular dynamic simulations have been un- 
dertaken in this work in order : i) to generate the amor- 
phous and ii) to simulate the annealing step and the re- 
crystallisation of an amorphous/crystalline silicon inter- 
face. In order to generate the disordered part, a melting 
process simulated by a thermal annealing at high temper- 
ature with a low quenching rate [14] is considered. The 
procedure has been repeated for the three clusters built 
along directions [100], [110] and [111] and containing re- 
spectively 4928, 5120 and 5184 atoms. The cell size pa- 
rameters (8 x 8 x 10) is a trade-off between reasonable 
computational costs and a possible influence of the peri- 
odic conditions on regrowth rate. Fig [T] illustrates the var- 
ious initial configurations. In order to simulate the layer 
re-growth, the Molecular Dynamics (MD) technique can 
also be applied to estimate the various atomic displace- 
ments at a given temperature in the framework of a Nose- 
Hoover thermostat [15] . The two most used empirical po- 
tentials Tersoff [H] and Stillinger- Weber (SW) [17] are 
under study in this work. The regrowth velocity deter- 
mination is based on the evolution of the c-Si/a-Si or /-Si 
interface. This location is estimated by using a plane struc- 
ture factor as a function of the depth [18] as shown in Fig. 
[2j The temperature range for annealing is a function of the 



empirical potential and is given in Table [T] Only lateral 
periodic boundary conditions are considered leaving the 
possibility to the stack to expand or contract. In order to 
perform a uniform doping in the amorphous part, a spe- 
cific approach has been developed. The amorphous part 
is sliced in various trenches of different thicknesses. Boron 
atoms are introduced in the amorphous lattice in a sub- 
stitutional way. After the dopants introduction, the equi- 
libration of the cluster is performed by a short annealing 
step leading to a final dopants distribution located in the 
amorphous part. Different concentration levels of boron 
are considered from the lowest one (0.2%) up to (5.4-6%). 
In the case of Molecular Dynamics with doped boron sys- 
tems, only the Tersoff potential has been used. It is often 
claimed that no empirical potential describes properly the 
subtle Si-B and B-B interactions. In this work, the ap- 
proach of Perez-Martin et al. [19] is followed where boron 
atoms are assumed to adopt a similar behaviour in silicon 
like carbon atoms. 



2.2 Recrystallisation velocity and disordered part 
properties 

In a previous work, it has been shown that for the simula- 
tion of silicon solid phase epitaxy the empirical potential 
must be taken carefully into account with the high sensi- 
tivity to describe or preserve an amorphous structure dur- 
ing annealing [20 . However, there is a debate on the influ- 
ence of the initial properties of the disordered part and the 
recrystallisation velocity obtained by molecular dynamics. 
In this work, the quality of the initial interface has been 
enhanced in comparison to previous studies where only 
a Wooten Winer Weaire's amorphous cluster on a crys- 
talline box was considered. The two main usual methods 
for generating a virtual amorphous cluster (Monte- Carlo 
bond switching [21] or by melting-quench) are confronted. 
The melting temperature used to generate a liquid part 
and next the quenching rate is described in Table [U As 
shown in Fig. [3] and also in table [2j large differences for 
the different properties (from radial and angular distri- 
bution to the mean coordination number distribution) in 
the disordered part can be observed. Fig. [4] reports the 
recrystallisation velocity of the a-Si or /-Si/c-Si interface 
for the two main potentials Tersoff and Stillinger- Weber 
(SW). For the recrystallisation velocity obtained using a 
Tersoff potential, Arrhenius kinetics have been observed 
for the two different amorphous in relative close agree- 
ment characterised with the same activation energy of ~ 
2.9 eV. The situation is completely different in the case 
of a simulated annealing using Stillinger- Weber. Despite 
the fact that in the first case, an initial a-Si/c-Si structure 
whereas in the second case a completely different /-Si/c-Si 
interface is under consideration, small differences can be 
observed. Clear liquid phase epitaxy characteristics with 
velocity re-growth of about a few ~ m/s in a wide tem- 
perature range limited by the melting of the crystalline 
part is observed slightly above a temperature of 1675K. 
These results strengthen our previous conclusions about 
the major influence of the empirical potential used for the 
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estimation of the recrystallisation velocity and on the vari- 
ability of empirical potentials to preserve an amorphous 
structure upon annealing [20] . 



3 Results and discussion 

3.1 Orientation effects and solid phase epitaxy 

3.1.1 Regrowth mechanisms 

The influence of the interface orientation effect on the re- 
growth mechanism has been investigated. The interface 
migration process is studied for the different crystalline 
orientations [100], [110] and [111]. Fig. [5] presents the a- 
Si/c-Si for an interface migration at 1700K during ther- 
mal annealing. For [100], the re-growth front is governed 
by first a diffusion process in the disordered part, then by 
a partial ordering at the interface and finally by a one- 
by-one atom incorporation. The propagation of the crys- 
talline front proceeds in a quasi-planar a-Si/c-Si interface 
with some partial reordering which could extend to one or 
two planes. 

For the orientation [110] a relative different mechanism 
has been observed in Fig. [6] Two silicon atoms have to 
be incorporated at the recrystallisation front and then to 
bind with each other. In order to get a second atom and to 
complete the crystalline plane, the recrystallisation is ob- 
served to propagate also vertically in another nucleation 
site. The re-growth front predicted by simulation is made 
of nanofacets leading to the formation of a sawtooth in- 
terface of about 3-4 atomic planes. For [110], the interface 
is more diffuse and characterised by a more important 
roughness as a partial ordering extend to several atomic 
planes. 

For the case [111], the recrystallisation mechanisms are 
observed to be even more complex (Fig. [7]) in comparison 
with [100] and [110]. The migration from a-Si to c-Si oc- 
curs by the rearrangements of complete or partial bilayers. 
For example, in the case of Fig. [7] incomplete bilayers form 
on the right and left part, while the central part is slowly 
rearranging (10 ns and 20 ns) letting suddenly appear sev- 
eral crystalline layers (25 ns). A buried defect at the top 
surface remains. The interface migration is made of suc- 
cessive plateaus whose height correspond to the distance 
(or a multiple of the distance) between the bilayers (111). 
The formation of one or more twins in planes (111) has 
been observed in nearly all the temperatures simulated. 
Other defects are observed in about 50 % of the simula- 
tions. The [111] interface recrystallisation interface is also 
completely different and mostly based on a bilayer process. 
For [111], the silicon atoms tends to organise partially in 
bilayer with a significant time necessary for the atomic 
arrangements within the the bilayer as three atoms of the 
bilayer are involved in the atomic process for the creation 
of a six- fold rings. Again, in the case [111], the recrystalli- 
sation process is not strictly planar with a large interplay 
between in and out of plane nucleation. 



3.1.2 Orientation dependent recrystallisation velocity 

One main open question concerns the influence of the re- 
growth mechanisms related to orientation effects at the 
a-Si/c-Si interface. Fig. [8] presents an example of the inter- 
face velocity extraction in the intrinsic case for the three 
orientations under study. This corresponds to the case 
where it was possible to define an a-Si/c-Si interface using 
a planar structure factor. This is not possible for in-plane 
and extended defects which generate too much disorder 
to allow a proper interface definition. All the systems con- 
sidered are free from defects except for systems [111] con- 
taining only twin parallel defects. In this particular case, 
the orientation [111] is estimated to be kept since it is only 
a default in the stacking plane sequence. The results as a 
function of the temperature and for the three orientations 
are given in Fig. [8] in the temperature range that could be 
considered with a reasonable computing time. It can be 
observed that while being very close in the orientations 
[100] and [110], the velocity is lower along the direction 
[111] by a factor 1.43. For the orientation [100], the acti- 
vation energy can be extracted using the classical theory 
of thermally activated growth [27]. Considering the fact 
that the interface velocity reported in the Arrhenius plot 
is superimposed or either parallel, there was no reason to 
consider a different activation energy for the cases [110] 
and [111]. 

3.2 Discussion 

The theoretical picture given by molecular dynamics sim- 
ulations stresses the importance of the initial interface ori- 
entation on the growth mechanisms. It has been observed 
that solid phase epitaxy mechanisms are strongly depen- 
dent on the interface orientation with a quasi-planar mi- 
gration for the [100] interface. This propagation is compat- 
ible with most analytical modelling based on a sharp in- 
terface migration. For [110] and [111], the situation is rela- 
tively different with the formation respectively of nanofacets 
for [110] and a bilayer mechanism leading to the formation 
and the migration of a more complex and diffuse inter- 
face associated with the generation of defects. Even if the 
physical reality described by molecular dynamics is more 
complex, some connections could be made with the phe- 
nomenological model of Drosd and Washbrun [22] which 
described the initial recrystallisation steps [23]. The re- 
sults on the interface recrystallisation velocity are more 
complicated to confront with experimental results. One 
main issue in our study is to define an interface migration 
velocity in the case where the a-Si/c-Si interface is more 
diffuse and characterised by the presence of defects. In this 
paper, only free-defect re-growth is considered. It corre- 
sponds to an ideal physical picture which is not observed 
experimentally [22] . The resulting activation energy of the 
process for the various orientation is 2.87 eV in accordance 
with the experimental range of [2.4-2.9]eV [22124125] . This 
activation energy cannot be associated with a simple ac- 
tive process like self-diffusion or defects migration. This 
energy seems rather associated with the energy barrier at 
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the interface as suggested by Posselt et al [26] . Since some 
bond breaking and exchange between the amorphous and 
the silicon phase and exchange are often observed in our 
simulation in order to allow a proper definition and propa- 
gation of the a-Si/c-Si interface. Some links could be made 
with the picture given by molecular dynamics and the ap- 
proach of Lu et al [27] which considered the dangling 
bond breaking, moving and trapping mechanism as the 
main process responsible for the growth rate limitation 
[8]. The prefactor for orientation [111] (1.2 xlO 7 m/s) is 
reduced by a factor 1.43 compared to [100] and [110] sug- 
gesting that the atomic incorporation mechanism based on 
three atoms at interface [111] slows down the overall in- 
terface migration compared to the other orientation. Only 
a small part of the velocity reduction can be explained by 
our molecular dynamics simulation. Since experimentally 
more interface recrystallisation velocity reduction is ob- 
served for [110] and [111], the role and the influence of 
defects on the nucleation velocity should be investigated 
more in depth at the theoretical level. A remaining open 
question that should be addressed concerns the influence 
of the cell size particularly in a large presence of defects. 
It seems to be the most plausible hypothesis to explain 
the reduction velocity effects. 



3.3 Regrowth of a doped a-Si:B/c-Si interface 

The possible influence of boron atoms in the silicon amor- 
phous part has been investigated in the case of interface 
[100]. The implantation step has not been Fig. [9] presents 
an example of initial doped a-Si (doped)/c-Si interface 
[100] and illustrates the final state where the amorphous 
part is completely recrystallised with most of the boron 
atoms incorporated in substitutional positions of the crys- 
talline matrix. This result is in agreement with a previous 
study on the annealing of damage generated by octodecab- 
orane implantation [28] . The results concerning the inter- 
face recrystallisation velocity are first detailed just below 
and then discussed in further subsection. 



3.3.1 Interface recrystallisation velocity 

Fig. [TUlreports the evolution of the interface position at 
2OOOK0 for different boron concentrations. In this case, it 
is clearly shown that the presence of boron is increasing 
the velocity of the interface migration. Fig. [TT1 summarises 
all the computer experiments performed in the intrinsic 
case and for 2 % and 6 % boron concentration configura- 
tions. Basically, for most temperatures, it can be observed 
that there is no significant difference between the intrinsic 
case and the 2 % case. The calibration by an Arrhenius 
law for the last case gives the same activation energy and a 
small difference for the prefactor (6.2 %). For a very large 
boron concentration (6 %), a small increase in the solid 
phase epitaxy velocity is observed leading to an increase 
in the prefactor by a ratio 1.34. 



3.3.2 Discussion 

In this work, a complete different behaviour has been ob- 
served compared to similar previous molecular dynamics 
simulations performed by Mattoni et al. [29 . Using molec- 
ular dynamics in the framework of the Stillinger- Weber 
empirical potential, an increase of the re-growth velocity 
at low boron concentrations (below 1%) has been observed 
whereas for larger concentrations (1% up to 10%) the ve- 
locity is reduced. On the experimental side [24], it has been 
observed that below 0.06 %, the boron concentration has 
almost no influence on the recrystallisation velocity. How- 
ever, around 0.2 %, a significant increase is observed up 
to a factor 10 at a low temperature annealing (~ 600°C). 
Our results obtained using Tersoff are much more in agree- 
ment with the experimental picture and the observed mi- 
nor influence of boron at low concentration. Probably, in 
the case of Stillinger- Weber simulation, the mechanism is 
more compatible with a framework of liquid phase epitaxy 
with the high tendency of SW based a-Si for melting even 
at low temperature [20] . Compared to the experimental 
results, it would not be so surprising that the influence of 
dopant is largely underestimated by our empirical molec- 
ular dynamics simulation since the Fermi level variations 
in our system are not taken into account theoretically. 
Generalised Fermi-Level Shifting (GFLS) model and the 
influence of charged defects at the interface are often as- 
sumed empirically with some success to estimate the en- 
hancement of the normalised re-growth rate defined by the 
ratio of the doped recrystallisation velocity in the doped 
with the intrinsic one VspE [3Q] : 



VSPE 
Vint 



Vint 



1 + g exp 



(Ef ~ Eg) 

kT 



g exp 



(Efi - E d ) 
kT 



(2) 



1 which would correspond to a rescaled temperature of about 
1000°C 



where Ef and are respectively the Fermi level and 
the defect level in the band gap whereas g is the degen- 
eracy factor. As pointed out by Johnson et al. [30] more 
work is needed in this direction in order to confirm the 
validity of the GFLS approach. First, a more refined de- 
scription level which incorporates electronic structure and 
Fermi shift such as tight-binding or ab-initio molecular dy- 
namics would probably be mandatory for such a purpose. 
Next, the ideal system under study (doped a-Si:B/c-Si in- 
terface) has probably to be improved in order to match the 
experimental picture since for example dopant implanta- 
tion is known to generate large amounts of defects (end of 
range defects, interstitials ,...) and surface roughness. 



4 Conclusions 

Empirical molecular dynamics simulation have been un- 
dertaken in order to investigate the silicon re-growth at 
the atomic level. Earlier results on the influence of the 
empirical potential and the large sensibility of the virtual 
silicon amorphous to the atomic force description have 
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been again illustrated. On the other hand, it has been ob- 
served that the interface migration velocity estimated is 
not so sensitive to the way a virtual amorphous silicon 
cluster is generated provided the amorphous character is 
preserved during annealing. The influence of orientation 
effects has been investigated at the atomic level showing a 
large influence of the surface orientation effects on the re- 
crystallisation mechanism. The physical picture described 
by molecular dynamics is more complex than the classical 
re-growth theory since the relation between the atomistic 
structure of the interface and the defects generation is 
clearly underlined. The concept of a sharp interface often 
used in the classical theory must be put in perspective 
as molecular dynamics suggests a more diffuse interface. 
More attention is probably needed to investigate more in 
depth the influence of defect on velocity re-growth. Con- 
cerning doping effects, a small increase in the recrystalli- 
sation velocity has been observed at large doping concen- 
tration. Further work is clearly needed in this direction 
to study if the influence of charged defects could be in- 
troduced empirically in molecular dynamics simulation in 
order to validate or not the generalised Fermi-level shifting 
model often assumed by the classical re-growth theory. 
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lar support on information technology aspects and to Fabrizio 
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Table 1. Main molecular dynamics parameters for disordered part generation. 



Potential 


Annealing Range 


Melting (K) 


Quenching 




(K) 




Rate 








(K/ps) 


Tersoff 


[1700-2500] 


3500 


1 x 10 11 


SW (Standard) 


[1000-1700] 


2200 


5 x 10 11 
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Table 2. Distribution of the mean coordination number (CN) 
the different initial disordered interfaces at 300K. 



Conditions 


CN 


2fold 


3fold 


4fold 


5fold 


6fold 


a-Si, WWW 


3.63 


7.6 % 


19.7% 


70.8 % 


0.9% 


0% 


a-Si, Tersoff 


4.17 


0.7 % 


9.9% 


63.4 % 


24 % 


1.9% 


/-Si, SW 


2.92 


24. % 


38.5% 


25.67% 


3.4 


0.16 % 
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Fig. 1. Stick models of the different a-Si/c-Si interface for 
the three orientations [100], [110], [111]. 
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Fig. 2. Example of in-plane structure factor indicating 
the location of the a-Si/c-Si interface. The horizontal line 
at Sq = 0.5 highlights the location of the a-Si/c-Si. 
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Fig. 3. Radial and angular distribution properties for the 
different potentials and the amorphous reference gener- 
ated by the Wooten Weiner Weaire approach. 



C. Krzeminski, E. Lampin: Some insight from molecular dynamics simulation on amorphous silicon re- growth 11 



E 





















-•- WWW+Tersoff 

www+sw 

— •— Tersoff+Tersoff 










— - sw+sw 















0.5 0.6 07 0^8 0^9 1 



1000/T(K" 1 ) 

Fig. 4. Evolution of the [100] interface re-growth veloc- 
ity as a function of the annealing temperature for the two 
main types : Wooten Weiner Weaire (WWW) approach 
and melting/quench and for the two main potentials Ter- 
soff and Stillinger- Weber (SW). The variation remains rel- 
atively modest compared to the difference in atomic clus- 
ter size (2000 atoms for WWW) and this work (~ 5000 
atoms) and those with the initial interface structure. 
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Fig. 5. Evolution of the a-Si/c-Si [100] interface position 
during solid phase epitaxy at 1700K. The evolution is re- 
ported for respectively 2.5, 5, 10, 15 and 20 ns illustrating 
the recrystallisation process where a quasi-planar inter- 
face propagates by a one atom by one atom incorporation 
mechanism. 
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Fig. 6. Evolution of the a-Si/c-Si [110] interface position 
during solid phase epitaxy at 1700K. The evolution is re- 
ported for respectively 5, 10, 15, 20 and 25 ns illustrating 
the recrystallisation mechanism where the formation and 
the dissolution of nanofacets is observed. 
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Fig. 7. Evolution of the a-Si/c-Si [111] interface position 
during solid phase epitaxy at 1700K. The evolution is re- 
ported for respectively 5, 10, 15, and 20 ns illustrating the 
bilayer recrystallisation mechanism and defect formation 
during re-growth. 
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Fig. 8. Evolution of the interface recrystallisation veloc- 
ity along orientation [100], [110] and [111] for the Tersoff 
potential. The orientation influence concerns mostly the 
prefactor since the activation energy of the regrowth pro- 
cess is observed not to be orientation dependent. 
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Fig. 9. An example of an initial a-Si:B/c-Si [100] inter- 
face for a 6 % boron concentration (top panel). The next 
figure shows the configuration after annealing where the 
dopant has been incorporated during re-growth. The seg- 
regation performed by the interface and the diffusion in 
the amorphous part can be clearly observed. 
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Fig. 10. Example of interface migration for various 
dopant concentrations at 2000K as a function of the an- 
nealing time. 
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Fig. 11. Recrystallisation velocity for solid phase epitaxy 
as a function of temperature and boron concentration in 
the initial amorphous phase. 



